Mesoscopic fluctuations of the Density of States and Conductivity in the middle of 

the band of Disordered Lattices 
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The mesoscopic fluctuations of tlie Density of electronic States (DoS) and of the conductivity of 
two- and three- dimensional lattices with randomly distributed substitutional impurities are studied. 
Correlations of the levels lying above (or below) the Fermi surface, in addition to the correlations 
of the levels lying on opposite sides of the Fermi surface, take place at half filling due to nesting. 
The Bragg reflections mediate to increase static fiuctuations of the conductivity in the middle of the 
band which change the distribution function of the conductivity at half- filling. 
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Recent studies of electronic level statistics in disor- 
dered systems have shown [Q-Q that existence of repul- 
sions between the levels in metallic phase results in a 
realization of the Wigner- Dyson statistics, ||]. Energy 
levels in a sufficiently dopped d > 2 dimensional electron 
gas become uncorrelated in an insulator phase and the 
level distribution obeys to Poisson statistics, How- 
ever, overlapping of one-particle states with different en- 
ergies leads to level correlations which change their dis- 
tribution. 

Our recent studies of weak localization effects in two- 
dimensional (2D) square lattice and in three- dimen- 
sional (3D) cubic lattice with substitutional impurities 
have revealed that Bragg reflections (BR) due to com- 
mensurability of the electronic wavelength. A, and a lat- 
tice spacing, a, strongly change the localization picture, 
pO|-p^ . The DoS vanishes on the Fermi surface for non- 
interacting electrons in a 2D lattice and it acquires a 
small dip on the Fermi surface of 3D simple cubic lat- 
tice with approaching half-filling, ||l^,|ll[. Nevertheless, 
electron-electron (e-e) interactions give a positive quan- 
tum correction to the DoS, which compensates the 
Altshuler-Aronov's negative logarithmic corrections to 
the DoS of 2D systems, 14|. Therefore it is interesting 



to clarify how static fluctuations of physical parameters 
of 2D disordered lattice with nested Fermi surface are 
changed by BR as the half-filling is approached. 

In this paper we consider the effect of BR on the level 
statistics and on conductivity fluctuations of 2D and 3-D 
disordered lattices with a half- filled band. A particular 
characteristics of level spectra is the two-level correlation 
function. 
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2, 3) lattice calculated in the Born approximation: 
Po2 = i^^lyit ln(e_Fmm{To, j^j}), and poz = const- 

^2t3%a.3 Vkl - ( at |e| « t) with t and Tq being the 
tunneling integral for nearest-neighbor sites and the re- 
laxation time for clastic impurity scattering, respectively. 

By using the formula for the DoS, p(e) = 
2^1 / ^{^^(r, r'; e) - Gfl(r, r'; e)}, which relates p(e) to 
the retarded (G^) and advanced (Ga) Green's functions, 
i?(e, e') can be expressed as. 



dr / dr' 



where (...) 
tions. Pod 



means averaging over impurity realiza- 
is the DoS of the d-dimensional {d = 



2nvpod 

{{GA{r, r; e)G^(r', r'; e')> + (Gfl(r, r; e)Gfl(r', r'; e')) - 
-{Gr{v, r; e)GA{r', r'; e')> - (GA(r, r; e)GR{r\ r'; e')) 
-4Re{GR{r, r; e))i?e(Gfl(r', r'; e'))}, (2) 

where v is the 'volume' and s is the factor of spin 
degeneracy. Far from half-filling the correlators RA 
and AR in Eq.(^) give only contributions to the two- 
level correlation function, However, existence of 
the electron-hole symmetry for nested Fermi surfaces 
gives rise to considerable contributions of the RR and 
AA correlators to R{e,e') in Eq.(^). The Fermi sur- 
face of a d-dimensional lattice with the energy spectrum 
of e(p) — ^X^iLifl ~ cos(pia)] becomes nested at half- 
filling, when ep — dt, which permits an electron-hole 
symmetry, e{p + Q) — ep — ~[e(p) — ef], with respect 
to the nesting vectors Q = {±7r/a, 7r/a} for 2D and 
Q = {zLtt / a, ztir/ a, IT /a} for 313 lattices. New singular 
impurity blocks take place at half-filling with particle- 
hole symmetry, which are referred to as the 7r-Diffuson 
(Dtt) and the 7r-Cooperon (G^r). The 7r-Diffuson (tt- 
Coopcron) has a diffusion pole at large oc Q momenta 
differences (total momenta) and small total energies of 
the electron and the hole (of two electrons), |10 1^. 

As it is known, an interference between the self- inter- 
secting trajectories leads to the weak localization of an 
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electronic wave function, (see, |13|). The electron passing 
the loop (e.g. in Fig. la), which is formed on the trajec- 
tory due to multiple scattering on impurities with small 
tilted angles, in clock- and counterclockwise, reduces its 
transmission probability. However, for a nested Fermi 
surface each act of impurity scattering is accompanied 
by BR with large (^ tt) scattering angle (see Fig. lb), 
which strongly changes the weak localization picture of 
free electron gas. Worthwile to notice that the electron- 
hole symmetry effects in the 2D case strongly differ from 
that in the one-dimensional (1-D) case, where BR act as a 
destructive factor of localization and result in big effects 
(like the Dyson singularity in the ID DoS , |lq ]) only in 
the absence of forward scattering (see, e.g. [[l6[ ) reversing 
backward scattering to forward one. 





FIG. 1. a) Self- intersecting trajectory due to multiple scat- 
tering on impurities with small tilted angles in i — th act of 
scattering; (b)The same trajectory that is drawn in (a) with 
the exception that each scattering is accompanied by BR, re- 
sulting in angles Qi -l-vr.In the both figures the magnitudes of 
velocity vectors are chosen to be the same. 

The expression for the two particle impurity block 
C^(q, e + t') in the particle-particle channel due to Umk- 
lapp scatterings with the particles' energies e and e' is 
given as, 



^TTPodTo ^ 



(«') 



where the phenomenological parameter 7 is introduced 
to signify an inelastic processes rate. The vr-Diffuson 
-D7r(q, e + e') is also expressed by Eq.(^) with exception 
that q will now be the momenta difference of a particle 
and a hole with accuracy of the nesting vector Q. Notice 
that the 'normal' Cooperon and Diffuson blocks depend 
on the difference of the energies e and e' instead of sum 
in Eq.d). 



New diagrams (see. Fig. 2) appear at half-filling due 
to BR which give contributions to R{e, e') in addition to 
that coming from normal scattering in a diffusive system. 
The study of the diffusive regime assumes that the linear 
length L of the system is larger than the elastic mean free 
path I and smaller than the localization length, the latter 
of which is exponentially large. By considering the energy 
scales as the Thouless energy, Ec — hD / with the dif- 
fusion coefficient D = of d dimensional system, and 
the average level spacing, A = 

that ^ = hpodDL'^-^ = 



it is possible to see 
— is a dimension- 
less conductance. Since g > 1 for the metallic case, the 
diffusive system can be characterized by the condition 
A ^ Ec <ti fi/To- By summing up contributions of the 
diagrams in Fig. 2 the correlator i?(e, e') is expressed as. 



T^{-i\e- e'\ + Dq^ +-fy 



[(l-*r,|e + e'|+7To)2 + |(q02-l]2. 
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q" = Z]a=i9a and = {2Tr /aNa)na with 
(— 7Va/2) < Ua < {Na/2), and (3 is the Dyson index clas- 
sifying the orthogonal, unitary and symplectic ensembles 
with P = 1, 2, and 4, respectively, Far from half- 
filling where BR are suppressed, only the first term in 
the bracket of Eq. M) contributes to i?, B . 
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FIG. 2. First order additional corrections to the DoS cor- 
relator R{e,e') in the metallic regime due to BR from a) 
TT-Diffuson and b) 7r-Cooperon blocks. 

In the ergodic regime, |e ± e'| <C Ec, only the q = 
term needs to be retained in the summation over q in 
Eq.(i. So, 
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which shows that two levels on the opposite sides of the 
Fermi surface with energy difference |e— e'| > 7 repel each 
other and they attract at energy difference of |e — e'| < 7. 
On the other hand, two levels with energies |e -I- e'| > 7 
on the same side of the Fermi surface attract each other 
and they repel for the energies |e -t- e'| < 7. Furthermore, 
attraction of two levels on the Fermi surface (e = e' = 0) 
weakens with approaching half-filling and the correlator 



2 



i?(0, 0) reachs its 3/4 value at half-filling. Notice that far 
from half-filling the levels lying only on the opposite sides 
of the Fermi surface interact with each other. Additional 
interaction of the levels on the same side of the Fermi 
surface appears due to BR at half-filling. 

The correlator of two levels centered at e,, and and 
averaged in an energy interval of E < W, where W = 2dt 
is the band width, can be obtained by integrating i?(e, e') 
given by Eq.(^ over the energy interval E: 
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£2-(i?/2 + *7)2][e^2 
+ e',f-{E + z^)^] 



(i?/2 + Z7)2][e^^-(i?/2 + Z7) 



iE/2 
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Re In 



E^ 



x{9{-eoe',)--9{eoe',)}; 



\eolK\>E/2. 
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Far from half-filling the second contribution in the brack- 
ets of Eqs.(^ and vanishes. This case corresponds to 
the continuum model, However, the fact that interac- 
tions of levels lying on opposite sides of the Fermi surface 
do give contribution to {5p^^{E)5p^i [E)) has not been 
taken into account in |^. As it is seen from Eqs.(^)-(0) 
the two-level correlation function strongly depends on 
center of energy strip E even a correlation is considered 
in the same energy interval when to = e'o- A logarithmi- 
cal energy dependence of the variance {[5pe^{E)Y) takes 
place for a strip centered around the Fermi level: 



{[5p,^=,{E)f) = 
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E/2<j< E, 
7 < E/2, (8) 



where / is the parameter characterizing the BR: / = 1/4 
at half-filling and / = far from commensurate points. 
According to Eq. (||) the Dyson repulsion of levels for en- 
ergies E /2 < -f < E turns to attraction of levels for large 
energy distances E > 2j. 

For the diffusive limit, when E ^ Ec, summing over 
q in Eq. can be replaced by integration. As a result 
we get the following expressions for the DoS variance 
{[5p..=,{E)Y) : 



([(5pe„=o(i?)]')*:/ = 



(%/2^1)(l-/). 



3/2_ 



d^2. (9) 



In the d = 

{[Spe^ME)?) 



2 case linear contributions in E to 
in Eq. (p|) are completely cancelled and the 



fluctuations are not as strong as in 31? systems. This 
seems to be connected with the localized character of 
levels in 2D systems. 
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FIG. 3. The diagrams which contribute to the conductivity 
variance due to BR. The diagrams denoted by primes differ 
from the presented ones through the direction of the electron 
lines. (e),(/) are symmetric to (e),(f) with respect to the 
single impurity line; and (g), (h) are obtained from (g),(h) 
by interchanging the straight and dashed lines under single 
impurity lines. Diagrams, similar to {e)-{h') exist also in the 
Cooper channel which are produced from (d, d'). 

Far from half-filling when the Fermi surface is approx- 
imately spheric the variance of fluctuations in static con- 
ductivity is designated by the diagrams given, e.g. in Fig. 
4 of Rcf . 1^ . These diagrams have been calculated also in 
Jl7| , p^ |3| for continuum and isotropic systems. Other con- 
tributions to the conductivity variance exist in the lattice 
model under consideration at the commensurate points 
due to BR which come from the diagrams given in Figj^. 
The total contributions to the conductance fluctuations 
due to Normal and Umklapp scatterings on impurities 
are: 
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(10) 

Here we followed the notations in where the contri- 
butions from the diffusion coefficient and the DoS fluc- 
tuations to the conductance variance (SGafiSGj^^) were 
denoted by the temperature dependent coefficients G^, 
and Gp, respectively. The expressions for G|, and G^ 
can be presented as: 

^2 K 



G' 



D 



q ' 



and 



r2 - _ 



TT ' J 2T'' ^2T' ^ 



(11) 



(12) 



where f{x) ~ 3L££th^j_^ T^he second term in the brackets 

r—^ ' Sinn a: „ 

in Eq. ( |ll| ) comes from the diagrams given in Fig.ra due 
to TT-scatterings. However BR give no contribution to 
Gp. According to Thouless picture , only one-electron 
states lying in an interval of Ef, centered on the Fermi 
level give contributions to the conductivity. Contribu- 
tions to Gp from the levels correlated on the same side 
of the Fermi level seem to cancel each other. 

At small temperatures T <^ Ec the coefficients Gjj and 
Gp do not depend on temperature: 



Gi 



{s^/(i){e^/TT^nfbd and 



Gn 



(1 



ml (13) 



where hd is a constant, which depends on the system 
dimension. In the case when T Ec the values of G^ 
and Gp strongly differ from each other and depend on 
temperature: 



Gi 



and, 



Gl=-{l + f)Gp; d = 3 



due to overlapping states of different energy in diffusive 
systems transform the distribution function from Poisson 
function to Wigner- Dyson one. Umklapp scatterings on 
impurities give additional contributions to the off- diago- 
nal matrix elements. Therefore, the change in the distri- 
bution function due to Umklapp scattering is reasonable. 
Two-level correlations are sensitive to whether the levels 
attract or repel each other and to the relative position 
of these levels, either on the same side or on the oppo- 
site sides of the Fermi surface. However, the conductance 
fluctuations seem not to be sensitive to the character of 
level interactions and the level positions; both attraction 
and repulsion give similar contributions to {SGapSG-yfj_) . 
This fact seems to be the reason why Eqs.(^)- (^ for 
{[Sp.^iE)]^) and Eqs.(|l3|),(|l|),(|l5|) contain the factor / 
coming from BR in different way. 

Recently fabricated Geo-based novel field-effect de- 
vices allow one to control the band filling by changing 
the gate potential. Half-filling is reached for 3 electrons 
doping per Geo molecule, @. The possibility of doping 
the fuUerites by substitutional impurities, while preserv- 
ing the periodicity of the Bravais lattice, will allow the 
observation of the commensurability effects on the meso- 
scopic fluctuations at half-filling in these devices. 
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where ad is some coefficient, g]. As can be seen from 
Eqs.(0 and (H) the main contribution to the conduc- 
tance variance comes from the fluctuations of the diffu- 
sion coefficient, which are intensified at half-filling. 
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